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1 Introduction 



In this paper we study the global structure of the cut locus (set of points reached optimally by more than 
one geodesic) for the simplest sub-Riemannian structures on three dimensional simple Lie groups (i.e. SU{2)^ 
6'0(3), and SL{2)) namely, the left- invariant sub-Riemannian structure induced by their Cartan decomposition 
and by the Killing form. 

Let G be a simple real Lie group of matrices with associated Lie algebra L and Killing form Kil(-, •). Let 
L = k0p be its Cartan decomposition with the usual commutation relations [k, k] C k, [p, p] C k, [k, p] C p. 
If L is non compact we also require k to be the maximal compact subalgebra of L. The most natural left- 
invariant sub-Riemannian structure that one can define on G is the one in which the distribution is generated 
by left translations of p and the sub-Riemannian metric < •, • > at the identity is generated by a scalar 
multiple of the Killing form restricted to p. The scalar must be chosen positive or negative in such a way that 
the scalar product is positive definite. We call G, endowed with such a sub-Riemannian structure, a k p 
sub-Riemannian manifold. 

k p sub-Riemannian manifolds have very special features: there are no strict abnormal minimizers and 
the Hamiltonian system given by the Pontryagin Maximum Principle is integrable in terms of elementary 
functions (products of exponentials). More precisely, if we write the distribution at a point ^ G G as A(^) = 
^p, we have the following expression for geodesies parametrized by arclength, starting at time zero from go 

g{t)=goe^^'^^-^'e-^'\ (1) 

where G k, G p, and we have < Ap^Ap >= 1. Thanks to left-invariance, with no loss of generality we 
can always assume go to be the identity and we will do so all along the paper. 

In all three-dimensional cases (i.e. SU(2)^ S0{3) and SL{2))^ p has dimension 2, while k has dimension 
1. Writing p = span{pi,p2} where {^1,^2} is an orthonormal frame for the sub-Riemannian structure (i.e. 
< Pi^Pj ^ij) k = span{/c}, we can write Ap = cos{0)pi + sin(^)j92 and A^ = ck with G M/27r, c G M. 
The map associating to the triple (^, c, t) the final point of the corresponding geodesic starting from the identity, 
is called the exponential map: 

5'^ X R X R+ ^ G 
* {0, c, t) ^ Exp(^, c, t) = e(^'«+^-)^e-^^^ 

For three dimensional k p sub-Riemannian manifolds, the local structure of the sub-Riemannian spheres, 
cut loci and conjugate loci starting from the identity has been described by Agrachev (unpublished) and, due 
to cylindrical symmetry of the Killing form in the p subspace, it is very similar to the one of the Heisenberg 
group. Indeed, locally, the cut locus coincides with the first conjugate locus (i.e. the set where local optimality 
is lost) and it is made by two connected one-dimensional manifolds adjacent to the identity and transversal to 
the distribution, see Figure [H 

Cut locus and 
first Conjugate locus 



Sub-Riemannian sphere 



ing point 

Figure 1: Local structure of sub-Riemannian spheres, cut and conjugate loci for 3-dim k p sub-Riemannian 
manifolds. 

However the global structure of the cut locus was still unknown. Indeed, to our knowledge, no global 
structure of the cut locus is known in sub-Riemannian geometry apart from the one of the Heisenberg group. 
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The main result of our paper is the fohowing: 
Theorem 1 Let Ki^ be the cut locus starting from the identity. We have the following: 

• for SU{2), i^id is a maximal circle without one point (the identity). 

• for S0{?>), Ki^ is a stratified set made by two manifolds glued in one point. The first manifold is WF'^ , 
the second manifold is a maximal circle without one point (the identity). 

• for SL{2), i^id is a stratified set made by two manifolds glued in one point. The first manifold is M? , the 
second manifold is a circle without one point (the identity). 

For ah cases the one dimensional strata contains the cut locus appearing in the local analysis. 




with antipodal identification 



Figure 2: The cut loci for the k p sub-Riemannian manifolds 5/7(2), 6'0(3), SL(2). 



Notice that k p sub-Riemannian manifold SU (2) has the structure of CR manifold and it is a tight 
structure [3116]. 

Once the cut locus is computed, one can provide the expression of the sub-Riemannian distance from the 
identity. The following theorem gives the sub-Riemannian distance for SU{2). The proof, given in Section 
I5.1.H can be adapted to get similar results in the cases of 6'0(3) and SL{2). 



Theorem 2 Let g 

Id is 



a P 
-p a 



eSU{2), i.e. a,^GC, 



1 . Its sub-Riemannian distance from 



where arg {a) G [0, 27r] and '0(a) = t the unique solution of 

( 



) (27r - arg(a)) if p = 

if /?^0 ' 



(2) 



-f +arctan^^ 



tan(^2m^^ =arg(a) 



t€ 0, 



(3) 



This theorem and its analogs for SO (3) and SL{2) are useful to give estimates for the fundamental solutions 
of the hypoelliptic heat equation induced by the sub-Riemannian structure (O [121 [El [IS])- Moreover this 
Theorem can be seen as the answer, in the case of SU{2), to the question (formulated in [Mj) about the 
possibility of inverting the matrix equation ([T]), i.e., for every matrix g G SU{2), find a matrix A = -\- Ap, 
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with < Ap^Ap >= 1, solution to the equation g = gQe^^k+Ap)t^-Akt ^ If ^ 7^ then this equation has one and 
only one solution, otherwise it has more than one solution (indeed infinitely many, see Sections [3] and [5]) . 

Then we study the most natural sub-Riemannian structures on the lens spaces q) induced by the one on 
SU (2). The lens space q) (with p, q coprime integers, p, ^ 7^ 0) is the quotient of SU (2) by the equivalence 
relation 

{-k %)~{-k *)if3-€,«,„o.o,„„i.,.c„.ha.(-) = (^ :,)(;;). 

The lens spaces are three-dimensional manifolds, but they are neither Lie groups nor homogeneous spaces of 
SU{2), except for the case L(2, 1) S0{?>) . 

In the case of lens spaces we get that the cut locus is much more complicated with respect to those on 
SU{2) and SL{2). It is still a stratified set, but in general with more strata. The precise description is given 
in Section [Ql 

Sub-Riemannian structures on the lens space L(4, 1) are particularly interesting for mechanical applications 
and for problems of geometry of vision on the two-dimensional sphere. Indeed, L(4, 1) PTS^^ the bundle of 
directions of S'^ . These applications are the subject of a forthcoming paper. 

The structure of the paper is the following: in Section [2] we recall the definition of sub-Riemannian manifold, 
we state the Pontryagin Maximum Principle (that is a first order necessary condition for optimality for problems 
of calculus of variations with non-holonomic constraints) and we define the cut and conjugate loci. Then 
we define k p sub-Riemannian manifolds. In Section [3] we define k p sub-Riemannian structures on 
SU{2)^ S0{3)^ SL{2) and compute the corresponding geodesies and conjugate loci. In Section [H we give sub- 
Riemannian structures on lens spaces as quotients of the k p sub-Riemannian structure on SU{2). The core 
of the paper is Section O where we compute the cut loci and the sub-Riemannian distance. The general idea 
is the following: we first identify the prolongation of the cut locus arising locally, then we compute the part of 
the cut locus due to the symmetries of the problem and finally we show that there is no other cut point. 

2 Basic Definitions 

2.1 Sub-Riemannian manifold 

A (n, m)-sub-Riemannian manifold is a triple (M, A,g), where 

• M is a connected smooth manifold of dimension n; 

• A is a Lie bracket generating smooth distribution of constant rank m < n, i.e. A is a smooth map that 
associates to g G M a m-dim subspace A{q) of T^M, and V g G M we have 

span{[/i, [. . . [fk-ijk] . . ]]{q) \ h e Vec(M) and Mp) G A{p) y p e M} = T,M. (4) 

Here Vec{M) denotes the set of smooth vector fields on M. 

• gg is a Riemannian metric on A{q), that is smooth as function of q. 

The Lie bracket generating condition (jH) is also known as Hormander condition. 

A Lipschitz continuous curve 7 : [0,T] ^ M is said to be horizontal if j{t) G A(7(t)) for almost every 
t G [0,T]. Given an horizontal curve 7 : [0,T] M, the length of j is 

^(7) = ^ yjg-,it)m),i{t))dt. (5) 

The distance induced by the sub-Riemannian structure on M is the function 

d{qo, qi) = inf{/(7) I 7(0) = ^0, l{T) = gi, 7 horizontal}. (6) 

The hypothesis of connectedness of M and the Lie bracket generating assumption for the distribution 
guarantee the finiteness and the continuity of <i(-, •) with respect to the topology of M (Chow's Theorem, see 
for instance [3J). 
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The function <i(-, •) is called the Carnot-Charateodory distance and gives to M the structure of metric space 
(see UW)' 

It is a standard fact that /(7) is invariant under reparameterization of the curve 7. Moreover, if an admissible 
curve 7 minimizes the so-called energy functional 

E{l)= f g^(t)(7W,7W)^t. 

with T fixed (and fixed initial and final point), then v = Y^g^(^)(7(t), 7(t)) is constant and 7 is also a minimizer 
of /(•). On the other side a minimizer 7 of /(•) such that v is constant is a minimizer of E{-) with T = l{'y)/v. 

A geodesic for the sub-Riemannian manifold is a curve 7 : [0, T] ^ M such that for every sufficiently small 
interval [^1,^2] C [0,T], 7|[^^ ^ minimizer of E{-). A geodesic for which g^(^)(7(t), 7(t)) is (constantly) 
equal to one is said to be parameterized by arclength. 

Locally, the pair (A, g) can be given by assigning a set of m smooth vector fields that are orthonormal for 
g, i.e. 

A{q) = Span{F,{q), . . . , F^(g)}, g,(F,(g), Fj{q)) = S,j. (7) 

When (A, g) can be defined as in ([7j) by m vector fields defined globally, we say that the sub-Riemannian 
manifold is trivializable. 

Given a (n, m)- trivializable sub-Riemannian manifold, the problem of finding a curve minimizing the energy 
between two fixed points go, ^1 G M is naturally formulated as the optimal control problem 

q = Y^u^F,{q), ii, gR, / ^ ^/,'(t) ^ min, g(0) = go, q(T) = q^. (8) 

It is a standard fact that this optimal control problem is equivalent to the minimum time problem with controls 
ui,. . . ,Um satisfying ul^ Vu^<l. 

When the manifold is analytic and the orthonormal frame can be assigned through m analytic vector fields, 
we say that the sub-Riemannian manifold is analytic. 

In this paper we are concerned with sub-Riemannian manifolds that are trivializable and analytic since they 
are given in terms of left- invariant vector fields on Lie groups. 

2.2 First order necessary conditions, Cut locus, Conjugate locus 

Consider a trivializable (n, m)-sub-Riemannian manifold. Solutions to the optimal control problem ([8]) are 
computed via the Pontryagin Maximum Principle (PMP for short, see for instance [31 [HI [2T1 [24]) that is a 
first order necessary condition for optimality and generalizes the WeierstraB conditions of Calculus of Varia- 
tions. For each optimal curve, the PMP provides a lift to the cotangent bundle that is a solution to a suitable 
pseudo-Hamiltonian system. 



Theorem 3 (Pontryagin Maximum Principle for the problem (|8])) Let M he a n- dimensional smooth 
manifold and consider the minimization problem in the class of Lipschitz continuous curves, where Fi, i = 
1, . . . , m are smooth vector fields on M and the final time T is fixed. Consider the map ^ : T*M x M x — > M 
defined by 

m m 

^{q,\pQ,u) := < X,^UiFi{q) > ^pQ^uj(t). 

If the curve q{.) : [0,T] M corresponding to the control u{.) : [0, T] is optimal then there exist a never 

vanishing Lipschitz continuous covector A(.) : t e [0,T] X{t) G T*(^^-^M and a constant pq < such that, for 
a.e. t e [0,T]; 

i) g(t) = — (g(t),A(t),po,^t)), 
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ii) A(t) = - — (g(t),A(t),po,^iW). 

iii) ^(g(t),A(t),po,^t)) = 0. 

Remark 1 A curve q{.) : [0, T] ^ M satisfying the PMP is said to be an extremal In general, an extremal may 
correspond to more than one pair (A(.),po)- If an extremal satisfies the PMP with po 7^ 0, then it is called a 
normal extremals. If it satisfies the PMP with po = it is called an abnormal extremal. An extremal can be 
both normal and abnormal. For normal extremals one can normalize po = —1/2. 

If an extremal satisfies the PMP only with po = 0, then it is called a strict abnormal extremal. If a strict 
abnormal extremal is optimal, then it is called a strict abnormal minimizer. For a deep analysis of abnormal 
extremals in sub-Riemannian geometry, see (S] [15] . 

It is well known that all normal extremals are geodesies (see for instance [3j). Moreover if there are no strict 
abnormal minimizers then all geodesies are normal extremals for some fixed final time T. This will be always 
the case in this paper: indeed we are concerned with sub-Riemannian manifolds of dimension 3, defined by a 
pair of vector fields Fi and F2 such that for all q G M, Span{Fi{q)^ F2{q)^ [Fi{q)^ F2{q)]} = TqM^ i.e. the so 
called 3-D contact case., for which there are no abnormal extremals (even non strict). 

In this case from iii) one gets Uiit) =< A(t), F^(t) >, z = 1 . . . , m and the PMP becomes much simpler: a 
curve q{.) is a geodesic if and only if it is the projection on M of a solution {\(t)^q{t)) for the Hamiltonian 
system on T*M corresponding to: 

-J m 

satisfying i^(A(0), g(0)) ^ 0. 

Remark 2 Notice that H is constant along any given solution of the Hamiltonian system. Moreover, = ^ if 
and only if the geodesic is parameterized by arclength. In the following, for simplicity of notation, we assume 
that all geodesies are defined for t G [0, +00). 

Fix qo e M. For every Aq G T^^M satisfying 

H{Xo,qo) = l/2 (9) 

and every t > define the exponential map Exp(Ao, t) as the projection on M of the solution, evaluated at time 
t, of the Hamiltonian system associated with with initial condition A(0) = Aq and q{0) = qo- Notice that 
condition (O defines a hypercylinder A^^ S"^'^ x R^-^ in T*^M. 

Definition 1 The conjugate locus from qo is the set Cq^ of critical values of the map 

ExD- ^oXM+ - M 

^""P- {\o,t) ^ Exp(Ao,t). 

For every Aq G A^q^ let t{Xo) be the n-th positive time, if it exists, for which the map (Ao,t) Exp(Ao,t) is 
singular at (Ao,t(Ao)). The n-th conjugate locus from qo C^^ is the set {Exp(Ao, t(Ao)) | t{\o) exists}. 
The cut locus from qo is the set Kq^ of points reached optimally by more than one geodesic, i.e., the set 



Kq^ = < G M I 3 Ai, A2 G Ago, Ai 7^ A2, t G such that 



g = Exp(Ai, t) = Exp(A2, t), and 
Exp(Ai, •), Exp(A2, •) optimal in [0, 



Remark 3 It is a standard fact that for every Aq satisfying ([9|), the set T(Ao) = {t > \ Exp(A,t) is singular 
at (Ao,F)} is a discrete set (see for instance [3]). 

Remark 4 Let (M, A, g) be a sub-Riemannian manifold. Fix qo ^ M and assume: i) each point of M is reached 
by an optimal geodesic starting from qo; ii) there are no abnormal minimizers. The following facts are well 
known (a proof in the 3-D contact case can be found in [4J). 
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• the first conjugate locus C^^ is the set of points where the geodesies starting from lose local optimality; 

• if is a geodesic starting from go and i is the first positive time such that q{t) e Kq^ U C^^, then q{.) 
loses optimality in f, i.e. it is optimal in [0, t\ and not optimal in [0, t] for any t > t; 

• if a geodesic q{.) starting from qo loses optimality at f > 0, then q{t) G Kq^ U C^^; 

As a consequence, when the first conjugate locus is included in the cut locus (as in our cases, see Section [5]), 
the cut locus is the set of points where the geodesies lose optimality. 

Remark 5 It is well known that, while in Riemannian geometry Kq^ is never adjacent to go, in sub-Riemannian 
geometry this is always the case. See [2|. 

2.3 k©p sub-Riemannian manifolds 

For the sake of simplicity in the exposition, all over the paper, when we deal about Lie groups and Lie algebras, 
we always consider that they are groups and algebras of matrices. 

Let L be a simple Lie algebra and Kil(X, Y) = Tr{adx o ady) its Killing form. Recall that the Killing form 
defines a non-degenerate pseudo scalar product on L. In the following we recall what we mean by a Cart an 
decomposition of L. 

Definition 2 A Cartan decomposition of a simple Lie algebra L is any decomposition of the form: 



Definition 3 Let G be a simple Lie group with Lie algebra L. Let L = k p be a Cartan decomposition o/L. 
In the case in which G is noncompact assume that k is the maximal compact subalgebra o/L. 

On G, consider the distribution A{g) = gp endowed with the Riemannian metric gg{vi^ V2) =< g~^vi^ g~^V2 > 
where < , >:= a Kil|p( , ) and a < (resp. a > 0) if G is compact (resp. non compact). 

In this case we say that (G, A, g) i5 a k p sub-Riemannian manifold. 

The constant a is clearly not relevant. It is chosen just to obtain good normalizations. 

Remark 6 In the compact (resp. noncompact) case the fact that g is positive definite on A is guaranteed by 
the requirement a < (resp. by the requirements a > and k maximal compact subalgebra). 

Let {Xj} be an orthonormal frame for the subspace p C L, with respect to the metric defined in Definition 
[3l Then the problem of finding the minimal energy between the identity and a point ^1 G G in fixed time T 
becomes the left- invariant optimal control problem 



This problem admits a solution, see for instance Chapter 5 of [13]. 

For k p sub-Riemannian manifolds, one can prove that strict abnormal extremals are never optimal, 
since the Goh condition (see [31) is never satisfied. Moreover, the Hamiltonian system given by the Pontrya- 
gin Maximum Principle is integrable and the explicit expression of geodesies starting from the identity and 
parameterized by arclength is 



where A/^ G k, G p and < Ap^Ap >= 1. This formula is known from long time in the community. It was 
used independently by Agrachev [1], Brockett [14J and Kupka (oral communication). The first complete proof 
was written by Jurdjevic in [22]. The proof that strict abnormal extremals are never optimal was first written 
in [10]. See also [lES]. 

Remark 7 In the 3-dimensional case, the Hamiltonian system given by the Pontryagin Maximum Principle is 
indeed integrable even if the cost is not built with the Killing form (bi- invariant), but is only left- invariant. For 
the case of SO {3) see 0. 



L = k p, where [k, k] C k, [p, p] C k, [k, p] C p. 



(10) 




(11) 
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3 SU{2)^ S0{?>)^ SL{2)^ their geodesies and their conjugate loci 

In this section we fix coordinates on SU{2), 6'0(3), SL{2)^ and we apply formula (pTj) in order to get the 
explicit expressions for geodesies and conjugate loci. 

3.1 The k © p problem on SU{2) 

The Lie group SU (2) is the group of unitary unimodular 2x2 complex matrices 

SU{2) = { ( ^ ) e Mat(2, C) | \a\' + = 1 

It is compact and simply connected. The Lie algebra of SU (2) is the algebra of antihermitian traceless 2x2 
complex matrices 

su{2) = 1^ -ia^ ^ Mat(2,C) | a G M,;^ G c| . 

A basis of su{2) is {^1,^2, k} where 

whose commutation relations are [^1,^2] = k [p2,k] = pi [k^pi] = p2. Recall that for su{n) we have 
Kil(X,y) = 2nTr(Xy), see [20l p. 186, 516]; thus for su{2) Kil(X, F ) = 4Tr(Xr) and, in particular, 
Kil{pi,Pj) = —2Sij. The choice of the subspaces 

k = span{/c} p = span {pi,p2} 

provides a Cartan decomposition for su{2). Moreover, {^1,^2} is a orthonormal frame for the inner product 
< •, • >= — ^Kil(-, •) restricted to p. 

Defining A (^) = gp and gc,('Ui, 'z;2) =< 9~^vi^ g~^V2 >, we have that (5'/7(2), A, g) is a k0p sub-Riemannian 
manifold. 

Remark 8 Observe that all the k p structures that one can define on SU(2) are equivalent. For instance, 
one could set k = spanjpi} and p = span{p2, k}. 

Recah that SU(2) ^ = | ^ ^ G | |ap + = l| via the isomorphism 

SU(2) 

a_ P \ fa 

-f3 a J ^ [p 

In the following we always write elements of SU (2) as pairs of complex numbers. 
3.1.1 Expression of geodesies 

We compute the explicit expression of geodesies using the formula (pTj) . Consider an initial covector A 
A(^, c) = cos(^)j9i + sin(^)p2 + c/c G Aid- The corresponding exponential map is 

Exp(l9,c,t) := Exp(A(l9,c),t) = e(^°«(^>i+«i"(^>2+c/c)^g-c/ct ^ 



csin(f)gy+?|) ^ ^^^^^^ cOs(vTT^f ) + Z (^ ^cos(f)gy+?|) _ ^.^^^^ cOsiVTT^I) 

(cos(f +^)+zsin(f + ^)) 
We have the following symmetry properties 
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cylindrical symmetry: 



• central symmetry: 



Exp((9, c, t) 



set 



Exp(l9, c, t). We have Exp(l9, -c, 



1 

e'^ 



Exp(0,c,t); 



K^-arg(/3))^ 



if p^o 

if /3 = 



3.1.2 Pictures of S'^ and 6'^ 

We recall a standard construction for representing 5^ in a two dimensional space and in a three dimensional 
one. For more details see e.g. [26]. Consider S'^ C and flatten it on the equator plane, pushing the 
northern hemisphere down and the southern hemisphere up, getting two superimposed disks D'^ joined along 



?4. 



it can 



their circular boundaries. The construction is drawn in Figure [3]- left. Similarly, consider C 
be viewed as two superimposed balls joined along their boundaries. In this case the boundaries are two spheres 
5^. A picture of is given in Figure [3]-right. 




Identified points 




Identified points 



Figure 3: Left: construction of the 2-dim picture of S'^. Right: the 3-dim picture of S^. 



3.1.3 The conjugate locus 

Recall that all the partial derivatives of Exp evaluated in (^, c, t) lie in TgSU (2) = g- su(2) with g = Exp(^, c, t). 



One can easily check that the three vectors g ^ • ^ff^| 
dependent (hence ^ is a conjugate point) if and only if 



^1 aExp 

^ l(^,c,t)' 



-1 ^Exp 
dt I 



(e,c,t) 



G su(2) are linearly 



2 sin (^\/l + c2 - Vl + cH cos ^v^l + c^^^ = 0. 



9 



The first term is if and only if^Ge^=|^Q^ ||a| = l|, while the second vanishes if and only if 
Vl + c^l = tan (Vl + c^l) 7 hence we have two series of conjugate times: 

• first series: t2n-i = JiJ^^.'^ to which correspond the conjugate loci Cl^~^ = \ Id; 

• second series: t2n = where {xi,X2,...} is the ordered set of the strictly positive solutions of 
X = tan(x), to which correspond the conjugate loci 



C2n 
Id 



' / csin( 



§^ (sin [i^) +^cos (^)) + cos(.„) (cos (^) - zsin (- 



c e . 



Remark 9 Notice that all the geodesies have a countable number of conjugate times. 

We present some images of conjugate loci (Figures H] and [5]) . For simplicity we present an image of their 
section with the plane Re {P) = 0. The complete picture can be recovered using the cylindrical symmetry. 



Figure 4: k p problem on SU{2): projection of the odd conjugate loci. 




Remark 10 Notice that the second conjugate locus is a 2-dimensional submanifold of SU{2), while the other 
even conjugate loci have self-intersections. 

3.2 The k © p problem on S0{3) 

The Lie group SO (3) is the group of special orthogonal 3x3 real matrices 

S0{3) = {g e Mat(3,R) | gg^ = ld,det{g) = 1} . 

It is compact and its fundamental group is Z2. The Lie algebra of SO (3) is the algebra of skew-symmetric 
3x3 real matrices 



so{3) = 






—a 


b 


a 





—c 


-b 


c 






G Mat(3,R) 
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A basis of so{3) is {pi,P2, k} where 






pi = u u -i P2 



whose commutation relations are [^1,^2] = k [p2^k] = pi [k^pi] = p2- Recah that so{3) and su{2) are 
isomorphic as Lie algebras, while SU{2) is a double covering of S0{3). 

For 50(3) we have Kil(X, F) = Tr(XF) so, in particular, K\\{pi^pj) = —2Sij. The choice of the subspaces 

k = span{A:} p = span {^1,^2} 

gives a Cartan decomposition for so{3). Moreover, {^1,^2} is an orthonormal frame for the inner product 
< •, • >= — ^Kil(-, •) restricted to p. 

Defining A (^) = gp and gc,('Ui, 'z;2) =< g~^vi^ g~^V2 >, we have that (6'0(3),A,g) is a k0p sub-Riemannian 
manifold. As for SU{2)^ all the k p structures that one can define on S0{3) are equivalent. 

3.2.1 Expression of geodesies 

Consider an initial covector A = A(^, c) = cos(^)pi + sin(^)p2 + c/c G Am: using formula (pTj) . we have that the 
exponential map is 

Exp(l9, c, t) := Exp(A((9, c), t) = e(^°«(^)^^i+«i"(^)^^2+c/c)tg-c/ct ^ 

/ Ki cos(ct) + K2 cos(26' + ct) + K3C sin(ct) Ki sin{ct) + K2 sin(26' + ct) - K3C cos(ct) K4 cos{e) + K3 sin(6') \ 
_ / -Ki sin(ct) + K2 sin(26' + ct) + K3C cos(ct) Ki cos(ct) - K2 cos(26' + ct) + K3C sin(ct) -K3 cos(6') + K4 sin(6') \ 

I cosf \/l + c'^t)+c'^ i 

\ K4 cos{e + Ct) - K3 sin{e + ct) K3 cos{e + ct) + sin(6' + ct) ^ ^^^^ ^ y 

l+(l+2c2)cos(yT+^t) ^_ l-cos(yT+^t) sin(yT+^^) ^_ c(l-cos( VT+^t) ) 
With i^i = ^ 2(T+^^^ ^' ^2 = 2(1+^2^, ^3 = '^1+,. % = ^ tV^2 

The set of geodesies has symmetry properties similar to the SU (2) case. The conjugate locus can be obtained 
from the one of the SU{2) by the canonical projection SU{2) S0{3). As for SU{2), all the geodesies have 
a countable number of conjugate points. 

3.3 The k © p problem on SL{2) 

The Lie group SL(2) is the group of 2 x 2 real matrices with determinant 1 

SL{2) = {g e Mat(2,R) | dei{g) = 1}. 

It is a non-compact group and its fundamental group is Z. The Lie algebra of SL{2) is the algebra of traceless 
2x2 real matrices 

sl{2) = l( ^ ^ GMat(2,R)|. 



c 



A basis of sl{2) is {pi,P2, k} where 



1/1 





' 






' -1 


) 








. 1 



P' = 2yo -1 

whose commutation relations are [^1,^2] = —k [P2^k] = pi [k^pi] = p2. For sl{n) we have Kil(X, y) = 
2nTr(Xy'), see hence for sl{2) Kil(X, F) = 4Tr(Xy') and, in particular, Kil{pi,pj) = 2%. The choice of 
the subspaces 

k = span{/c} p = span {^1,^2} 

provides a Cartan decomposition for sl{2). For sl{2) the Cartan decomposition is unique, since k must be 
the maximal compact subalgebra. Moreover, {^1,^2} is a orthonormal frame for the inner product < •, • >= 
^Kil(-, •) restricted to p. 

Defining A(^) = gp dind gg{vi^V2) =< g~^vi^g~^V2 >, we have that (6'L(2),A,g) is a k0p sub-Riemannian 
manifold. 
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3.3.1 Expression of geodesies 

Consider an initial covector A = A(^, c) = cos{0)pi + sin(^)p2 + c/c G Aid: using formula ([TT]) . we have that the 
exponential map is 



Exp((9, c, t) 



Exp(A(i9,c),t) 



g(cos(0)pi+sin(0)p2+c/c)^g — cfet 



_ / cos (c|) + (cos (6' + c|) + csin (c|)) 

~ \^ -Ki sin (c|) + K2 (sin (6> + c|) + ccos (c|)) cos (c|)'+ ( 



Ki sin (c|) + K2 (sin (9 + c|) - ccos (c|)) \ 
cos (6> + c|) + csin (c|)) J 



with 



K2 



Cosh (vT 
cos I 



Sinh(VT^l) 



sin(yH2^|) 



[-1,1] 

c e (-00,-1) u (i,+oo) 

cG(-l,l) 
CG{-1,1} 

c e (-00,-1) u (i,+oo) 



3.3.2 A useful decomposition of SL(2) 

Proposition 1 For every g G SL(2) there exists a unique pair r G e^, s G eP such that g = rs. 



Take r 



cos(^) -sin(^) \ ^ k . 
sin(^) cos(^) ^"""^^ 



G SL{2). Notice that (0^,7) 7^ (0,0). We 



Proof. First, notice that = 6'0(2) and eP is the set of 2 x 2 symmetric matrices with determinant 1 and 
positive trace. 

a -\- S 13 — ^ 
/3 + 7 a — 5 

have to prove that exists a unique G M/27r such that s = r~^g is symmetric with positive trace. By direct 
computation one gets that s is symmetric if and only if asin(6>) = 7cos(6>). For any (a, 7) G M^\(0,0) there 

cos(6>i) sin(^i) 
-sin((9i) cos(l9i) 

cos(^2) sin(6>2) 
-sin(6>2) cos(6>2) 

trace, either Oi or O2 provide Tr(5) > 0. ■ 
Topologically 5'^ and eP :^ R^, hence SL{2) x R^. In the following, we represent SL(2) as the 

set R^ X [0, 1] with the identification rule (a, 6, 0) ~ (a, 6, 1). See Figure [H 



exist two solutions of this equation ^1,6^2 ^ R/27r with ^2 = ^1 + tt. Thus Tr 

-Tr I I ^^^'^^^^ ^^^l^^x I ^ ) . Observing that a symmetric matrix with determinant 1 has nonvanishing 




Identified planes 



Figure 6: A picture of SL{2). 



3.3.3 Symmetries in the SL{2) problem 
We have the following symmetry properties: 

• cylindrical symmetry: Exp((9,c,t) = e^o/cg^pi+yp2 ^^ere ( ^ ) = ( si^n(^) ~cos{f)^ ) ( ^0 ) 
(xo,yo,^o) are defined by Exp(0,c,t) = e^o^e^o^^+^o^^ 
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. central symmetry: Exp(^,-c,t) = e--o/c,.p.+,p. ^^^^^ ( ^ ) = ( si'n^^^ -co%l) ) ( ^ ) ^"""^ 
(xo, ^0, ^o) are defined by Exp(6>, c, t) = e^o^e^o^^+^o^^ 

3.3.4 The conjugate locus 

With similar arguments to those of Section [3.1. 3[ one checks that g = Exp(6>, c, t) is a conjugate point if and 
only if 



Sinh (vT^^i) (2Sinh (vT^^|) - tvT^^Cosh (vT^^|)) =0 c G (-1, 1) 

±^ = c = ±l 

sin(Vc2T7Ti) (2sin(Vc2T7T|) - tVc^^cos (Vc^^D) =0 cG (-oo, -1) U (1, oo) 

The first 2 equations have only the trivial solution t = 0. The third one gives two series of conjugate times: 

• first series: t2n-i = J^^i which correspond the conjugate loci C^^~^ = \ Id; 

• second series: t2n = where {xi,X2,...} is the ordered set of the strictly positive solutions of 
X = tan(x), to which correspond the conjugate loci 



cos{xn) cos{yn)-\-^^P==icos{e)-\-csin{yn)) cos{xn)sin{yn)-\-^^p^^{sin{e)-ccos{yn)) \ , f G 

Id 



2n ) I Vc2-i Vc^-i 

I I - cos{xn)s\n{yn) + ^^^=^{sin{e)+ccos{yn)) cos(x^) cos(2yri) + ^y=^(- cos(6')+csin(2yn)) / ^ G M/27r 



with = 



Remark 11 Notice that not all geodesies have conjugate points. Indeed, Exp(6>,c, •) has a conjugate point if 
and only if c G (—00, —1) U (1, +00). 

We present an image of the 2^^ conjugate locus (Figure [7j). For simplicity we present an image of its intersection 
with The complete picture can be recovered using the cylindrical symmetry. 




Figure 7: k p problem on SU{2): section of the 2^^ conjugate locus. 
Remark 12 Notice that all even conjugate loci have self-intersection. 

4 A sub-Riemannian structure on lens spaces 
4.1 Definition of L{p^q) 

Fix 2 coprime integers p^q G P^Q ¥^ ^- The lens space q) is defined as the quotient of SU (2) w.r.t. the 
identification rule 

I ~ f I if 3 G C p-th root of unity such that f ] = f ^ ) ( 

Lens spaces are 3-dimensional compact manifolds, but excepted L{2^1) :^ 6*0 (3), they are neither Lie groups 
nor homogeneous spaces of SU{2). The following topological equivalences hold: \/ p^q^k G Z, p^q coprime, 
g 7^ we have L{p, q) :^ L{p^ —q) ^ L{—p^ q) L{p^q -\- kp). Lens spaces have highly non-trivial topology, 
for details we refer to [25] . 

The following theorem permits to choose a representative of q) in SU{2). 
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Proposition 2 Consider the set 



eSU{2)\ Re (a) > 0, 



sin(^)' 



< 1 > C 5/7(2) and define 



dE^ = dEp n {Im (a) > 0}^ 9^^^ = dEp fi {Im (a) < 0}. Endow Ep with the equivalence relation ^ defined as 
follows: 



1. the relation is reflexive: 



G dE^ and ^ ^_ 



G dEp , we have 



2. moreover, given 

• either: Im(a+) = -Im(a") ^ anti /?+ = e^^'^i (5~ ; 



• or; Im (a+) = Im (a ) = {) and = e^^^ p (3 for some n e {1^ ... ^p}. 
The manifold Ep/^ is diffeomorphic to L{p^q). 



Proof. Take 



G 5'/7 (2) and let us look for uj p-th root of unity such that 



is equivalent to 



Re {coa) > and ^"^ ^"^""^ + \u;^P\^ < 1. 



^ ^ , e Ej). This condition 



(13) 



sm 



Recalling that \u;^l3\'^ = = 1 — |ap and that Im {coa) = \a\ sin(arg(a;a)) if a 7^ 0, equation ([13]) is equivalent 
to 



arg(cja) G 



TT TT 



or a = 0. 



Thus: 



if a 7^ 0, there exist at least one solution uoi of arg (coa) G 



TT TT 



(14) 



. Moreover, we have 2 distinct 



solutions Cc;i,cj2 if and only if arg (coia) = —f and arg (0020^) = f • ^^is case 



ujia 



\a\e p 



and 



\a\e p 



observe that 



uoia 



UJ2OL 

'I 



if a = 0, every co p-th root of unity satisfies 



^ ^ Ep] observe that for all the pairs uji , UJ2 we have 







Hence V 



G SU (2) we have a unique 



uoa 



L{p,q) = SU{2)/. 



G Ep/ i.e. the function 

uoa 



is bijective. 

Remark 13 A crucial observation for what follows is that the projection 



n : 



SU{2) ^ L{p,q) 
9 ^ [9] 
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is a local diffeomorphism. Moreover, H]^^ : Ep ^ L{p^ q)\ [dEp] is a diffeomorphism. In particular, Ep contains 
only 1 representative for each equivalence classes of q); i.e. if g^h e Ep and [g] = [h], then g = h. 

Remark 14 Proposition [2] provides a picture of L{p^q): recall that SU{2) is drawn as 2 balls in (see 
Section [3.1.2p . Hence Ep C SU{2) is drawn as a closed ellipsoid inside one of the 2 balls, via the map p : 
Bi{0) C 

. The picture of Ep is 



Ep 
a 



(Re(/3),Im(/?),Im(a)) 



(a;i,a;2,a;3) e Bi(0) 



sm 



< 1 



the one of Ep is Fp, see Figure [SI- left. The identification ^ induces the following identification on Fp-. given 
(x^,X2",X3") G dF^ = n {x3 > 0} and (xj",x^,X3 ) G dF~ = n {xs < 0}, they are identified when 



4 



cos(^) — sin(6>) 
sin(6>) cos(6>) 



with 6> = see Figure [H^right. 




Full ellipsis 




Identified points 




Figure 8: Left: F4. Right: the representation of L(4, 1), with some examples of the identification rule. 



Remark 15 Observe that the identification rule on Fp gives a 1-to-l identification between dFp fi {^3 > 0} and 
dFp n {xs < 0}, while there are in general more identified points on {x^ + = l} H {xs = 0}, see Figure [H 



Figure 9: L(4, 1): some examples of the identification rule on {xf -\- x^ = l} f) {xs = 0}. 
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4.2 Sub-Riemannian quotient structure on L{p, q) 

Proposition 3 The sub-Riemannian structure on SU{2) given in Section lSAl induces a 2-dim sub-Riemannian 
structure on L{p^q) = SU(2)/ via the quotient map 

SU{2) ^ L{p,q) 

X ^ [x] 

i.e. 

• the map _ 

A : [^] n* (A(/i)) C T^g]L{p, q) with h G [g] 

is a 2-dim smooth distribution on L{p^ q) that is Lie bracket generating; 

• ^[g](v^,w^) =< v^,w^ ^[g]'-^^ with h G [g], v,w G ThSU{2), U^{v) = v^, Il^{w) = is a 
smooth positive definite scalar product on A. 

Proof. The role of the map 11 and 11* |^ is ihustrated in the fohowing diagram 

TgSU{2)^^T[g]L{p,q) 

SU{2) ^L{p,q). 

The map E is a local diffeomorphism, thus n*|^ : TgSU{2) T^gjL{p^q) is a linear isomorphism, hence 

(A(^)) is a 2-dim subspace of T^g^L{p^q). 
The following statements: 

• the distribution A([^]) is well defined, i.e. V/ii,/i2 ^ [g] we have n*|^ {A{hi)) = n*|^ (A(/i2)) 

• the positive definite scalar product < v^^w^ >[^] is well defined, i.e. V/ii,/i2 ^ [^], vi^wi G Th^SU{2)^ V2^W2 G 
ThoSU(2) such that U^. (vi) = II*, (V2) and E*, (wi) = E*, (W2) we have < vi.wi >hi=< 

V2,W2 >h2 

are consequences of the following lemma: 
Lemma 1 Let /ii, /12 ^ [g] with /12 = ^ q ^ ^i* -^^^ 

P ^ P 

m \ ( Re(cc;^-i) -lm(cj^-i) \ / ni 



mi Im(cj^-i) Re(cj^-i) mi 

0/ \ o/\o 



is bijective. Moreover, it is an isometry w.r.t. the positive definite scalar product < , > and satisfies V 7^ G p 



dt\t^o dt\t^ 



Proof. Let = ( ^ ) G 6'/7(2) and 7? = (n, m, 0) G p. We have 



a cos [s/n^ + ^^|) - ^sin (Vn^ + m^f"^ 



1^ 6cos(Vn2 + m2|) +asin(Vn2 + m2|) 
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Take /ii,/i2 ^ [g] with h2 



u 
cc;^ 



hi and 771,772 ^ P with coordinates rji 



(77,2,7712,0). Consider the trajectories [hie 



trji 



and 



LV 



a cos 
6 cos 



6 sin 



sm 



(711, 7711, 0) and 7/2 

12; 



^2iA ni-\-imi 



[^2e 



/ 
LV 



ija cos 
co^b cos 

a cos 
6 cos 



^2i 
^2 2 



^2i 
h 2 



cj^^sin 
f cja sin 



t2 t 
^2 2 



t2 t 
^2 2 



A n2 + 



''2 2 



sm 



''2 2 



771 



2 t 
22. 



771: 



2 t 



712— im2 
n2+im2 



Thus 



le 



t77l 



[/i2e^^^] in the case 



cj^ -^(771 + imi) = 712 



77J + nil = ^2 + ^2 
77i — imi — uo^~^{n2 — 17712) 
77i + imi = uj^~^{n2 + 17712) 
17712. This equation is verified for 772 = (j^ivi)- 



22; v^^^fT^ J 

that is equivalent to 



Since 11 is a local diffeomorphism, V ^ G SU{2) 3 5(^) such that the map n*|^^ ^ : TB(g)SU{2) 

TB{[g])L{p^q) is a diffeomorphism, thus A is smooth. Lie bracket generating, and < v^^W:^ >[^] is smooth as a 
function of [g\. ■ 
Proposition [3] implies that the sub-Riemannian structures on SU{2) and L{p,q) defined above are locally 
isometric via the map 11. As a consequence, the geodesies of g), A,g) are the projection of geodesies of 
(S'/7(2), A, g). The conjugate locus for L{p^q) can be obtained from the one of SU{2) by the projection 11. 

Remark 16 One can check that the sub-Riemannian structure induced by SU{2) on L(2, 1) S0{3) is equiv- 
alent to the k p sub-Riemannian structure on SO (3) defined in Section [321 



5 Cut loci and distances 

In this section we prove the main theorems of the paper, i.e. we compute cut loci for SU (2), 50(3), lens spaces, 
SL{2) and we prove the formula (|3|) for the sub-Riemannian distance on SU{2). 

Recall that our problems satisfy the following assumptions: i) each point of M is reached by an optimal 
geodesic starting from Id, see Section 12. 3t ii) we are in the 3-D contact case, thus there are no abnormal 
minimizers. Hence Remark |4] applies. 

Proposition 4 Let T{0^ c) he the cut time for Exp(6>, c, •) (possibly +00 i/Exp(^, c, •) is optimal on [0, +00) ). 
Define V = {(i9, c, t) G Am x R+ | < t < T(i9, c)} and M' = M\ {Ki^ U Id). The function Exp,^ : V ^ M' 
is a diffeomorphism from V to M' . 

Proof. Let us first check that Exp(P) C M' . By contradiction, let Exp(6>,c,t) G M\M' , thus either t = 
or t = TiO^c) or Exp(^,c, •) is not optimal in [0,t], i.e. t > T{0^c). Contradiction. Let us verify that Exp|^ 
is injective: by contradiction, let Exp(6>i, ci, ti) = Exp(^2, C2, ^2) with (6>i,ci,ti) 7^ (6^2,^2,^2). If ti 7^ ^2, 
one of the two geodesies Exp(6>i, ci, •), Exp(6>2, C2, •) has already lost optimality, thus ti > T{6i,Ci)^ hence 
(Oi^Ci^ti) P, contradiction. If ti = ^2, we have that Exp(^i, ci, ti) is a cut point, hence ti > T(6>i,ci), 
contradiction. To verify that Exp|^ is surjective, take g G M' and observe that there is an optimal geodesic 
Exp(^, c, •) reaching it at time t < T{0^ c). But t = T(^, c) implies g G i^id, thus t < T{0^ c). 

The smoothness of Exp|^ and of its inverse follows from the fact that Exp is a local diffeomorphism outside 
the critical points (i.e. points where the differential of Exp is not of full rank) and that the critical points do 
not belong to V. Indeed, by contradiction, let (O^c^t) G P be a critical point, hence t is a conjugate time: 
it is either the first conjugate time, that coincide with the cut time (i.e. t = T{6^c)) or a greater one (i.e. 
t > T{6,c)). In both cases (O^c^t) V. Contradiction. ■ 
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5.1 The cut locus for SU{2) 

Theorem 1 The cut locus for the k © p problem on SU{2) is 

Ku = e^\U={e''' |ce(0,47r)}. 

Proof: Let ^Ge^\Id=|^ ^ ^ |aGC, |a| = l, a7^l| and Exp(^, c, •) the minimizing geodesic steering 

Id to g in time T. As a consequence of the cyhndrical symmetry, we have that Exp(?/^, c, T) = g V G ]R/27r, 
thus ^e^\Id) C Km. 

The core of the proof is to show that there are no cut points outside e^. Recah the expression of geodesies 

Exp{e,c,t) = ( ^ 

c.in(^)sin^^l) ^ ^^^^^^^^^^^—^^^^ ^ . ^ „os(^)WIT?^) _ ^.^^ ^ ^ ^^^^ , ^ 

.in(v^4) ^^^^^^ ^ + isin(f + 6)) 

By contradiction, assume that ^ G 5'/7(2)\e^ is reached by two different optimal trajectories Exp(^,c, •) and 
Exp(?/^, d, •) at time T. Observe that Exp(6>, c, y=^), Exp(?/;, ^^^) G C i^M, thus 

< T < min I . . ^ I . (16) 



(17) 



Observe that Exp(6>, c, T) = Exp('0,G^, T) imphes that \(3\ is equal in the two cases, i.e. 

sin(^I±2T) sin(^E^) 



From this equation it follows \c\ = \d\. Indeed equation (fT7|) is equivalent to 1.^ — = A.^ ^ . From 

2 2 

the fact that Vi+c^^ ^ ^/ i+d'^T ^ j^q^ ^j^^^ ^j^^ function is injective for p G (0,7r), it follows that 

Thus we consider the two cases: 

• c = d e R: the cylindrical symmetry implies either that 6> = (so the 2 geodesies coincide) or ^ G e^. 
Contradiction. 

• c = —d G ]R\{0}: with no loss of generality we assume c > 0. Since by the central and cylindrical 
symmetries we have 

Exp(V^,-c,t) = ^ gi(^+^-^arg(/3))^ ^ whcrc Exp(l9,c,t) = ^ 
the equation Exp(6>, c, t) = Exp(?/^, — c, t) implies Im (a) = 0. Hence 



ct\ . /vTTc^A f ct\ I VTT~^t 



ccos — sm 



2/12/^ V 2 



yi^sin cos 11^^ . (18) 



The terms c, sin are non-zero because of equation (p!6|) and c < Vl + c^. 

Thus cos (^) = if and only if cos (^^^^^^^ = 0, that is impossible because < < ^-^^^^^ < 



TT. Hence we rewrite equation (fT8|) as — '^^^^^^ ^ = ^^"4^ state that (0,tan(0)), (^,tan(^)) , 

^ Vi+c^t ^ Vi+c^t ^^ g^^g 2 distinct points aligned on the graph of the function tan in [0,7r). It is 
impossible. Contradiction. 



(f) 
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The cut locus for the k p sub-Riemannian manifold SU{2) is given in Figure [TOl 



Identified points 




Identified points 



Figure 10: The cut locus for the k p sub-Riemannian manifold SU{2). 



5.1.1 The sub-Riemannian distance in SU{2) 

In this section we compute the sub-Riemannian distance on SU{2)^ i.e. we prove Theorem [2l 
^1 = 1 Q j G e : g is reached by a geodesic Exp(^,c, •) at time ^^^^2 



Let g 



Observe that Exp (O, c, 



= Exp l9, ± 



47r2 



COS 



TT^ — 1 q= z sm 



7r± 



. Thus the distance d{g^ld) is the smallest t > such that e 



iarg(a) 



whose solution is t = 2 y^arg (a) (27r — arg (a)) where arg (a) is chosen in [0, 27r]. 



Let ^ = 



/9 



G 6'/7(2)\e^. Applying proposition [H we have that Exp|J^(^) is well defined on V = 

|(^, c, t) G Aid X \ < t < -^==|. Thus the sub-Riemannian distance of g from the origin is d{g^ Id) = t 
where t is the third component of Exp^^(^), i.e. the unique solution t of Exp(6>, c, = g with {O^c^t) G P. 

Using the explicit form of Exp given in (I3.1.ip . one checks that the system | ^^^^^.l ^' ^ is equivalent to 



" 2 



h arctan 



: tan 



^^^"V 2 -)) =arg(a) 



cos(f + 6>) + i sin(f + 6>) = arg(;^) 

The third equation has no role for the computation of distance as a consequence of the cylindrical symmetry. 

Remark 17 The distance is a bounded function: this is due to its continuity and the compactness of SU{2). 
The farthest point starting from Id is —Id, whose distance is 27r. 
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Notice that V a, A, ^2 e C, \(3i\ = |/?2| we have d (^(^ j ,Idj = ^ (^(^ j 'J ^ ^ \\ ^1 J 'j' 
this is due to the cyhndrical and central symmetries. 

5.2 The cut locus for SO (3) and lens spaces 

In this section we compute the cut locus for lens spaces q). As a particular case, we get the cut locus for 
50(3)c^L(2,l). 

Theorem 4 The cut locus for the sub-Riemannian problem on L{p^q) defined in Section\^is a stratification 



with 



KIT = idE,] = 



[Id] 

k 



^ / X „ Im(a)'^ 
a,be C,Re(a) > 0, ^ ' 



T^IOC 



e 



\ m = 



sin (f ) 



Proof. Let us first prove the following lemma. 

Lemma 2 A geodesic in L{p^ q) steering [Id] to [g] in minimum time T admits a unique lift; 7o(-) in SU{2) 
starting from Id. 

Moreover, jo{t) = Exp(l9o, cq, t) W t e [0, T] for some Oq G M/27r, c G M. 

Proof. Take 7(-) as in the hypotheses. Since L{p^q) and SU{2) are locally diffeomorphic via 11, there is a 
unique lift 7o(-) in SU{2) starting from Id, i.e. 7o(0) = Id and [-fo{t)] = lit) V t G [0,T]. 

Let us prove that 7o(-) is an optimal trajectory reaching 7o(T). By contradiction, there exists a trajectory 
7i(-) such that 71 (ti) = 70 (T) with ti < T. Hence, its projection [7i(')] satisfies [71(^1)] = [g] with ti < T. 
Contradiction. 

Since 7o(-) is an optimal trajectory, it is a geodesic of SU{2)^ and there exist 6 G ]R/27r, c G M such that 
7o(t) =Exp(^o,co,t) G [0,T]. ■ 

Let us prove that C ^[id] • consider [g] G , a geodesic steering [Id] to [g] in minimum time T with 

unique lift Exp(6>o, cq, •)• % definition of K^^^^^, we have Exp(6>o, cq, T) G e^\Id C SU(2), i.e. Exp(6>o, cq, T) 
lies in the cut locus for the sub-Riemannian problem on SU{2). Thus there exists another optimal geodesic 
Exp(6>i, ci, •) defined in [0, T] such that Exp(6>i, ci, T) = Exp(6>o, cq, T) G [g]. Thus the geodesic [Exp(6>i, ci, •)] 
reaches [g] in minimum time. The geodesies in SU (2) are distinct in a neighborhood of Id, so their projections 
in a neighborhood of [Id] are distinct as well. 

Let us now prove that K^-^^ C i^[id]- consider [g] G K^-^^^ a geodesic steering [Id] to [g] in minimum 

time T with unique lift Exp(6>o, cq, cah Exp(6>o, cq, T) = ^ ^ ^ id]- If = we have [g] G K^^^^-^ or 

[g] = [Id], so assume P 0. Due to the cylindrical and central symmetries, we have Exp(6>o + V^, — co,T) = 

g2i(0o-arg(/3))+i^^ ^- Consider G M/27r solution of e2^(^o-^^s(^))+^^'^ = e^^^'p and i/j' G M/27r solution of 

g2i(^o-arg(/3))+#- _ e"^^^?. If Exp((9o , cq , T) G dE^, we have [Exp(6>o + -cq, T)] = [Exp(l9o, cq, T)] = [^^]; 
if Exp(^07Co,T) G we have similarly [Exp(^o + '^"7— co,T)] = [Exp(^0 7 cq, T)] = [^]. If cq 7^ 0, we have 

found two distinct trajectories reaching [g] in optimal time; if cq = 0, we have Exp(^07Co,T) G H 
thus at least one among ?/^+ and ?/^~ are not null, so at least one among Exp(^o + ^^^7 0? •) and Exp(^o + ^^"7 •) 
are distinct from Exp(^o, 0, •) in a neighborhood of Id, so are their projections in a neighborhood of [Id]. 

Finally, consider [g] G q)\ U K^-^^ U [Id]^ and assume by contradiction that there exist 2 distinct 

geodesies steering [Id] to [g] in minimum time T with distinct lifts Exp(6>o,co, Exp(6>i,ci, •). There are two 
possibilities: 
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• Exp(6>o, Co, T) = Exp(6>i, ci, T). In this case Exp(6>o, cq, T) lies in the cut locus for the sub-Riemannian 
problem on SU{2), hence [g] G K^^i^y Contradiction. 

• Exp(^07Co,T) 7^ Exp(^i, ci, T): since by hypothesis [g] ^[il^^ we have Exp(^0 7 cq, T), Exp(^i, ci, T) 
dEp. Recalling that [Exp((9o, cq, T)] = [Exp(6>o, cq, T)] and Exp(6>o, cq, T), Exp(6>i, ci, T) implies 
Exp(l9o,co,T) = Exp((9i,ci,T), due to remark [131 we have that Exp(l9i,Q,T) G _SU{2)\E^ for z = 
or z = 1. We assume with no loss of generality that Exp(^o,co,T) G SU{2)\Ep^ thus the geodesic 
Exp(l9o,co,t) with t G [0,T] connects Id G and Exp(l9o, cq, T) G SU{2)\E^- thus 3 t G (0, T) such 
that Exp(6>o, Co, t) G It implies 7o(?) = [Exp(6>o, cq, ?)] G ^[^j^p thus 7o(t) is no more optimal for 
t G [0,T]. Contradiction. 



Remark 18 Notice that i^T^^j^^^^ is a manifold (a circle without a point), while K^^^ is not in general. Indeed, it 
is an orbifold. It can be seen as S'^ C with the following identification: {x^ t) e S^n {X3 > 0} and 

(x7 ,x^) G 6*^ n ixs < 0| are identified when xf = —x^ and | I = I ^^^In''' I I ^5- 

^^ ^ '^^ Lo_j ^ ^ ^^^^^ sm^t/) cos^t/j J \ X2 

with l9 = In the case S0{3) L(2, 1), we have that K'^^ = RP^ (see Figure [U]- left), while in the other 
cases it is not locally euclidean: in fact, take a neighborhood of a point P on the equator and observe that it 
is topologically equivalent to a set of p half-planes with a common line as boundary. 

Next we give an idea of the topology of the cut locus for L(4, 1). Consider the space Ti made by the two 
intersecting strips {(a, 6,0) G | a, 6 G [-1,1]} and {(a, 0,6) G | a, 6 G [-1,1]} with the following identi- 
fication: (—1,6,0) ~ (1,0,6) and (—1,0,6) ~ (1,-6,0). The boundary of this set is topologically a circle S^. 
Consider now a 2-dimensional semisphere T2. The cut locus K^-^^ is topologically equivalent to the space given 
by glueing Ti and T2 along their boundaries S^. The cut locus i^Tfi^] is given by glueing K^-^^ with a circle 
along a point on T2 and next removing a point on (the starting point). See a picture of it in Figure [TTl-right. 



[Id] 





Identified boundaries 



with antipodal identification 




Figure 11: Left: The cut locus for the sub-Riemannian problem on S0{3). Right: the cut locus for the 
sub-Riemannian problem on L(4, 1). 
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5.3 The cut locus for SL{2) 

Theorem 5 The cut locus for the k © p problem on SL{2) is a stratification 
with 

Kir = e^-''eP = {geSL{2)\g = g^,T:tg<0}, 

Proof. Let us first prove that K^^^ C i^^id- Let ^ G \ Id: it is reached optimahy by a geodesic Exp(6>, c, •) 
at time T. Due to the cyhndrical symmetry, we have g = Exp('0, c, T) W i/j e ]R/27r, thus g G Ki^. 

Let us now prove that K^^^ C i^id- Let g = e^^^e^°^^+^°^2 ^ e^^^eP: it is reached optimahy by a geodesic 
Exp(6>, c, •) at time T. If + ?/q = we have g = e^^^ G K{^^, thus it is a cut point. If Xq^i/q ^ 0, 

due to the cyhndrical and central symmetry, we have Exp(6> + ip^—c^T) = ^-'^'^k^xpi+yp2 ^ith ^ ^ 



. Choose ip in such a way that ^ + is the angle on the plane of the 



cos(2(9 + V^) sin(2(9 + V^) \ ( xq 
sin(2(9 + V^) -cos(2(9 + V^) / V 

line passing through (0,0) and (xo^yo). In this way we have ( ) = ( )• Observing that e~^^^ = 

V ^ y V ^0 y 

we finally have that g = Exp(6>, c, T) = Exp(^ + i/j^ — c, T). Observe that c 7^ because Exp(6>, 0, •) G eP, thus 
the two geodesies Exp(6>, c, •), Exp(6> + — c, •) are distinct. 

We prove now that there is no cut point outside K^^^ U K{^^. By contradiction, let g G SL{2) \ 
(iCj^^^ U i^j^^ U Id) be reached by two optimal trajectories Exp(^, c, •) and Exp(?/^, •) at time T. Writing 

Exp(^,c,t)=f ^-{^'^'JM we have 



' Sinh(VT^l) 



:= V (^11 - ^22)^ + {g\2 + ^21)^ = 



CG (-1,1) 



t CG{-1,1} . (19) 

(^i) cG(-oo,-l)U(l,+^) 



sm 



The identity Exp(^, c, T) = Exp(?/^, T) implies r(c, T) = r((i, T), that implies = dP . Indeed, observe that 
in the 3 cases described by (fT9|) we have respectively r(c, t) > t, r(c, t) = t, r(c, t) < t, thus either c, (i G (—1,1) 
or c,(i G { — 1,1} or c, (i G (— 00, — 1) U (l,+oo). In each of the 3 case the identity r(c, T) = r{d,T) implies 
= (P. Indeed we have the following 3 cases 

case c, d G ( — 1,1): in this case the conclusion follows from the fact that ^^"^^^^ is injective for p G (0, +00); 

case c, dG { — 1,1}: straightforward; 

case c, d G (—00, —1) U (1, +00): Let us prove first that ^^^2"^"^ ^ (O^'^)- contradiction, assume ^^'^2^^ — 
TT. There exists t G (0,T] such that = 0, hence r{c^t) = 0, from which it follows Exp(6>,c, t) G e^. 

Hence either t < T (and Exp(^, c, •) is not optimal on [0, T], contradiction) or t = T (and g G i^i^^ U Id, 
contradiction). Similarly one prove that e (0,^). 

Now observe that r{c,T) = r{d,T) implies = Recalling that is injective 



for p G (0, tt) we have = hence = 

We have two cases: 
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• c = d e R. The identity ^ii(6>,c, T) = ^ii(?/;,c, T) implies either 6 = i/j (i.e. the geodesies coincide) or 
c e (-00, -1) U (1, +00), thus sin (Vc^ - if) = 0^ i-^- Exp(6>, c, T) G e^, hence either ^ = Id or ^ G i^^g^ 
Contradiction. 

• c = G M\{0}. Writing g = e^^e^o^^+^o^^ ^ave Exp(V;, -c,T) = e-^^e^^i+^^^ The identity 
Exp(^, c, T) = Exp(?/^, — c, T) and the uniqueness of the decomposition (|3.3.2p imply e^^ = ±Id. Thus g 
is symmetric, i.e. gi2{0^c^T) = ^2i(^,c, T). 

xr / ^ ^\ 1 • • • T tan(c^) Tanh ( \/l — ^ ) . i i i n 

It c G (—1,1) this equation implies — = ^1,^2^ • Choosing c > 0, observe that the nrst 

positive solution Ti of the equation ^^"(^^ ) = Tanh(vi c ^ ) gg^|.-g£gg ^ (71-^3:1^ Xhe other cases 

Vl — ^ 

c G { — 1, 1} and c G (—00, —1) U (1, +00) are treated similarly and lead to Ti G (tt, 3f ). 

Thus cos (c^) < 0, hence Tr (g) < 0. But Exp(^, c, Ti) symmetric and Tr (^) < implies Exp(^, c, Ti) G 
K^^^ ^ i.e. Ti is a cut time. Thus either T = Ti (meaning that Exp(6>,c, t) G K^^^) or T > Ti and 
Exp(^, c, •) is not optimal in [0,T]. Contradiction. 

■ 

We give a picture of the cut locus for the k p sub-Riemannian manifold SL{2) in Figure ([T2|) . 
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